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ABSTRACT — In the First Part of this paper [that was submitted for pub. in 1991 and 
appeared in print in Phys. Reports 214 (1992) 339] we critically review and analyse the 
main theoretical definitions and calculations of the sub-barrier tunnelling and reflection 



o 

00 . 

Q\ . times. Moreover, we propose a new definition of such durations, on the basis of a recent 

general formalism of ours (already tested for other types of quantum collisions) within 
conventional quantum mechanics. At last, we discuss some surprising results regarding 
the temporal evolution of the tunnelling processes: namely, the fact that QM predicts 
that tunnelling through opaque barriers takes place with Superluminal group-velocities. 
Aims of the Second Part [that appeared in print later, in J. de Physique-I 5 (1995) 1351] 
are: (i) presenting and analysing the results of various numerical calculations (based 
on our equations) on the penetration and return times < r Pen >, < r Rct >, during 
tunnelling inside a rectangular potential barrier, for various penetration depths xf, (ii) 
putting forth and discussing suitable definitions, besides of the mean values, also of the 
variances (or dispersions) Dtt and Dtr for the time durations of transmission and re- 
flection processes; (iii) mentioning, moreover, that our definition < tt > for the average 
transmission time results to constitute an improvement of the ordinary dwell-time r Dw 
formula. The numerical evaluations confirm that our approach implied, and implies, the 
existence of Superluminal tunnelling (that we call "Hartman effect"): an effect that in 
these days (due to the theoretical connections between tunnelling and evanescent-wave 
propagation) is receiving — at Cologne, Berkeley, Florence and Vienna — indirect, but 
quite interesting, experimental verifications. Eventually, we briefly analyze some other 
definitions of tunnelling times. A more detailed review of the same topics (in Italian, how- 



ever) has been "published" electronically, as LANL Archives # |cond-mat/9802126| . At 
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last, a brief review on the experimental data, that — in four different sectors of physics — 
seem to indicate the existence of Superluminal motions, appeared as an Appendix to the 
paper in LANL Archives # |physics/9712051| (to be published in Physica A, 1998). 



PACS nos.: 73.40. Gk ; 03.80. +r ; 03.65.Bz . 
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FIRST PART 



I-l - INTRODUCTION. 

During the last years, many attempts were devoted to theoretically defining and 
calculating the time spent by a sub-barrier-energy particle for tunnelling through a po- 
tential barrier; so that critical reviews were already in order. Actually, a first, valuable 
review article, by Hauge and Stoveneng, did recently appear (at the end of 1989) as ref.W. 
We deem it useful, however, a second review paper since: (i) on one hand, it appears 
convenient to deepen and extend the criticism to the existing approaches, and (ii) on the 
other hand, we can show — as a consequence of such a criticism — how a definite proposal 
can be put forth for the introduction of suitable, self-consistent, physically meaningful 
definitions of the tunnelling times: a question that was regarded still as an open problem 
in ref.M. 

The problem of defining the tunnelling times has a long history, since it arose in 
the fortieth and fiftieth' 2-4 ' , simultaneously with the problem of a general definition in 
quantum mechanics of the collision durations. Let us explicitly recall that the motion in- 
side a potential barrier is a quantum phenomenon without any classical or quasi-classical 
analogues, so that it lacks an easy intuitive representation; this can explain the birth of 
a variety of approaches to describe it, sometimes in contradiction with one another. 

The advent of high-speed experimental devices, based on tunnelling processes in 
semiconductors, did revive an interest in the whole question; ' 5,1 1 whose relevance has 
always been apparent in nuclear physics: For instance, tunnelling plays an essential role 
in the physics of nuclear fission and fusion. 

As already mentioned, aims of this paper are a comparative analysis of the various, 
known definitions for the tunnelling durations; and the introduction of a possible, rigorous 
way for theoretically evaluating the time spent by a particle inside a potential barrier: 
with the final presentation of some peculiarities of the tunnelling-process time evolution. 
Our starting point will be the formalism in refs.' 6 ' 8 '. 

First of all, in Sect. 1-2, we critically review the main definitions of the tunnelling 
times appeared up to now, to the best of our knowledge. Then, Sect. 1-3 is devoted to 
the application of our own formalism t 6 ' 8 ' to tunnelling. At last, Sect. 1-4 is devoted to the 
conclusions and to some proposals for further developments. 



1-2 - ABOUT THE EXISTING DEFINITIONS OF TUNNELLING TIMES. 
1-2.1 The simplest stationary picture of Tunnelling. 
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Let us confine ourselves to the ordinary case of particles moving only along the re- 
direction, and consider a time-independent barrier V(x) in the interval (0, a): see Fig.I-1, 
in which — for later convenience — a larger interval (x±, X2), containing the barrier region, 
is also indicated. We assume the stationary scattering problem to have been solved 
exactly for every kinetic energy E = h 2 k 2 /2m, where k is the wave-number and m the 
particle mass. The wave function ip(x] k) will have the general form: 



ip = ipj = ip in + ip R for x < 0, 

ip = tp u for < x < a, (1) 
■0 = 1pm = ip T for x > a; 



where: tp- m = e lkx ; = e %kx \ ipT = A T e lkx ; the lower indices R,T staying for 
"reflected" and "transmitted", respectively. In the simple case of a rectangular barrier 



[V(x) = Vo], it is ipn = ae KX + (3e KX [with k = y2m(V — E)/h], where the coefficients 
[amplitudes] A R , A T , a and (5 can be analytically calculated, and are known to be 

2k k + in 2k —k + m 

a = ~Dl {k 2 -K 2 )D_/D + + 2tkK ] ' Ul (k 2 - k 2 )D_/D + + 2ikn ex ^~ 2Ka ^ 

(2) 

A R = a + f3-l; A T = (ae Ka + (3e Ka ) e~ ka ; D± = 1 ± exp(-2/ca). 



The amplitudes Ar, At satisfy the probability conservation law 



|A R | 2 + |A T | 2 = 1. (3) 

The flux density 

does not depend on x. Before the barrier it equals the difference (1 — |Ar| 2 ) between 
the incoming and reflected wave flux-densities. It is less known that, inside the barrier, 
the fluxes for the separate components of ipu (exponentially decreasing and increasing, 
respectively: ae~ KX and f3e KX ) do vanish. Only their interference does provide the conser- 
vation of j. 
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1-2.2 Construction of the wave-packet 

For later use in the non-stationary description of actually moving wave-packets, let 
us consider a wave-packet constructed in terms of the solutions ip(x; k) of the station- 
ary Schroedinger equation: namely, by integrating over E from to oo with a weight- 
amplitude g[E — E) 



V(x, t)= J dE g(E - E) ip(x; k) exp(-iEt/h) (5) 

where we introduced the resolution of the time-evolution operator, with the normalization 
condition / dE \g\ 2 = 1, quantity E being the average kinetic energy. 

In the case of a Gaussian wave-packet it is convenient to pass from the energy to 
the impulse representation, by recalling that dE = (h 2 k/2m)dk; when the spread in E 
of g(E — E) is much smaller than E, one easily gets 

h 2 Tc h 2 

g(E -E)dE^ G{k -k)dk = ^— g{^- [k 2 - k 2 ]) dk, (6) 

with G = C exp[— b(k — k) 2 }. Of course, the (initial) wave-packet of the incoming waves 
will have a Gaussian shape also in the configuration space. 

By inserting in the integral of eq.(5) the reflected (V'r) o r transmitted (V't) wave, 
instead of the total wave ip, we obtain the final reflected or transmitted wave-packets, 
respectively, carrying a time-delay due to the interaction. Notice that one could expect 
a distortion in the wave-packet form due to the energy dependence of and A^; but 
it has been already shown, for a wide class of weight amplitudes, such a distortion to be 
negligible.^ Furthermore, we shall get rid also of the wave components with above-barrier 
energies by introducing the additional transformation 

g ( E - E) — g(E - E) Q(E - V ), (7) 

in order to avoid distortions of the sub-barrier penetration (tunnelling). 

For simplicity's sake, we shall in general address ourselves to quasi-monochromatic 
packets, for which the energy spread AE is so much smaller than E that it is possible to 
adopt the approximation 

\g(E — E)\ 2 ~ 8(E — ~E) (8) 
in all the final expressions for our quantities, when averaged over pdx or Jdt; where 
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p=|^(:r,t)| 2 ; J = Re[—^{x,t) ^ ] (9) 

are the probability density [for a particle to be located in the unitary space-interval 
centered at x] and the probability-density flux [for a particle to pass through position x 
during the unitary time-interval centered at t], respectively. 



1-2.3 The ordinary phase-times 

Following the usual procedure, introduced in' 2-4 ! , it is easy to get the ordinary phase- 
times for quasi-monochromatic packets, in the stationary-phase approximation 

r^(x h x f ; E) = - {x, - Xi ) + h (10) 
v ah 

and 

rl\x^E) = \{-2x l)+ n^q^ (11) 

where v = hk/m is the group-velocity; X{ G region I, and Xf G region III. Eqs.(10) and 
(11) refer to a transmitted [from the initial position x- x to the final position Xf\ particle 
and to a reflected [from the initial position X{ to the same position] particle, respectively; 
cf., e.g., ref.W. 

For a rectangular barrier with height Vq, the phase-times (10) and (11), when linearly 
extrapolated W to the barrier region (xi — 0; x? — a) would become 



4 h (0, a; E) = rl h (0, a; E) = — — [2Kak 2 (K 2 - k 2 ) + A; 4 Sinh(2/ta)] , (12) 

which, for na » 1, would simply yield 2/vk. In eqs.(12), it is D = 4K 2 k 2 + /c 4 Sinh 2 (/ta); 
and ko = 2mVo/h. In other words, ^ for sufficiently wide — i.e., opaque — (or high) bar- 
riers, eqs.(12) do not depend on the barrier width a, and the effective tunnelling-velocity 
a/r Ph may become arbitrarily large [Hartmann and Fletcher effect' 9 ' 10 ']. 

One of the main objections against extrapolations (12) is that they do not describe 
the actual asymptotic behaviour of the phase-times; since they disregard the fact that 
both the [magnitude of the] initial packet mean-position, and quantity Xf — a (where 
X{ is the transmitted packet mean-position) must be large with respect to the packet 
spatial extension [of the order of hv/AE], in order to avoid "interference" effects between 
physically quite different processes (i.e., between incident and reflected waves). 
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Therefore, it is not completely correct to attribute to the extrapolated phase-times 
the physical meaning of "times spent in the barrier region (= inside the barrier)". More- 
over, one cannot separate in and r^ h the self-interference delays from the time spent 
inside the barrier. 

Before going on, let us clarify the behaviour of the phase-times at the very top of 
the barrier, and check whether there is any continuity — there — between the values of 
the sub-barrier tunnelling time and those for the above-barrier case. Let us compare 
eqs.(12) with the following expression for the above-barrier transmission time: 



Ph , _ 2m -(k 2 - g 2 ) 2 tan(gq) + 4qak 2 (k 2 + q 2 ) / cos 2 (qa) 

r T (U, a, Z > V ) - Mq 4k2q2 + [{k2 + g2)tan(ga)]2 W 

which was obtained^*- 1 , by the stationary-phase method, for the case of a rectangular 
barrier. In such a case, it is ipu = r )t iqx + 5e~ iqx with q = ^2m[E — V )/h, and the 
coefficients 7 and 5 are analytically evaluable. By comparing eqs.(12) and (13) one gets 



P u . „ » mko? 2ma 

hm Tt (0, a; E > V ) = -77 , „ . — ► -— - , 

«-o T v ' ' 0J 6h{l + k 2 a 2 /4) a-00 3hk ' 



pv, / ^ . mko? 2ma 

\imT^(0,a;E <V ) = —- 7^777 — ► 7-77- 

ic-o T v u; 6^(l + A; 2 a 2 /4) «-oo 3^ 



(12') 



;i3') 



In other words, we find that that the two limits (12'), (13') do coincide with each other, 
and linearly depend on a for "opaque" barriers (provided that the condition no — > 
holds). Notice that such a result does not contradict the Hartmann and Fletcher effect, 
since the latter takes place only when kcl — > 00, while it is absent for finite values of na. 



(*) These calculations have been explicitly performed by V.S. Sergeyev. 
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1-2.4 The dwell time 

The total scattering time duration has been defined in' 11 ] as the probability for the 
particle to be localized in the interval between the initial (maybe, source) position and 
the final (maybe, detector) position, divided by the incident particle flux density; that 
is to say, as the time spent by a particle while travelling inside such space-interval: the 
so-called dwell time. In the chosen case of particles moving only along x, the dwell time 
is therefore defined as^ 12 ' 

1 r x t 

T Dw (x u x f ,k) = - dx \^(x;k)\ 2 . (14) 
For a rectangular barrier, the (dwell) time spent inside the barrier becomes: ^ 



.Dw/n i\ mk ro / 2 ,2\ , 2mV 



(0, a; k) = — — \2KaU 2 - k 2 ) + — ^Sinh(2Ka)] (15) 

which, for kcl » 1, would give Tik/nVo. The results (12) and (15) are in sharp contrast 
with each other with regard to the /c-dependence. Let us comment on this point. 

The stationary definition (14) for the dwell time, according to us, is not self-consistent 
from its very beginning, and appears to be in contradiction with its physical meaning. 
In fact, the time variable is firstly discarded (in passing from the time-dependent to the 
stationary Schroedinger equation), and later on it is re-introduced in an artificial, ad hoc 
way: namely, through the introduction of a localization-probability expressed in terms of 
time-independent wave functions, instead of actually moving wave-packets. 

Moreover, even the modified "dwell time approaches" with time-dependent wave 
functions! 13-15 ^, in which 



=Dw , , _ /\ df/* 1 dr i>(.r. t) 



x h x f ; fc) = , (16) 

J-oo dt Jm{Xi,t) 

do still contain formal time-averages, that are not actual averages over the phys- 
ical time (i.e., the time t(x) at which the considered particle passes through the 
position x). In fact those averages should be obtained — at least when the di- 
rection of flux J is time-independent — by integrating Jdt, and not pdtJ 6 ' 16 ^ In 
eq.(16) quantity J- in is defined as in eq.(9), just replacing ip(x;k) of eq.(5) by 
ip- m (x; k) = e %kx . At last, the "dwell-time approaches" are unable' 6 ' to define 
and study the time distributions for any kind of collision process. 
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1-2.5 The local Larmor times 

In' 17 ' a gedanken experimente was proposed for measuring the scattering duration as 
the ratio 9/u, where 9 is the angle by which the magnetic moment fi of the considered 
particle is rotated due to a small homogeneous magnetic field B = B z (directed along 
z) supposedly present in the scattering region, and uj = 2/iB /h is the Larmor precession 
frequency. For a magnetic field existing in the interval (xi, Xf), and for an incident particle 
(moving along x and) with spin | polarized along the x-direction (see Fig. 1-2), 9 results 
to be proportional to the average spin component < s y >: namely, 6 = — 2< s y > T /h, or 
9 = — 2< s y > R /h, for the transmitted or reflected waves, respectively. In this case, the 
Larmor times Ty^(x[, Xf ; k) and Ty R (xi, xf, k) become equal to the Phase Times, plus 
—however — terms which do oscillate as kx{ and kx{ vary.' 1 ' In the particular case of a 
rectangular barrier one gets 

r^(0, a; k) = r^(0, a; k) = r Dw (0, a; k). (17) 

In Baz's approach, as it was shown in refs.^, the expressions for the collision du- 
ration [e.g., eqs.(17)] are artificially distorted by the sharp boundaries attributed to the 
magnetic-field region; in other words, are influenced by the mathematical, rather than 
physical, assumptions. Actually, the oscillating terms do depend on the kind of boundary 
[for instance, smoothed] that one adopts. Moreover, both those oscillating terms vanish 
when Xi — > — oo and X{ — > oo, once one does average over the incident particle energy- 
spread; so that the final expressions do coincide with the Phase Times. 

It it also known' 1 ' that the mathematical behaviour assumed for the magnetic-field 
boundary does not influence only the spin components along y. In fact (see Fig. 1-2), the 
incident particle has finite probabilities of being spin-up or spin-down along the field- 
direction z. As pointed out in ref.' 12 ', the spin-up components will be preferentially 
transmitted (except when d|Ar| 2 /di? < 0): so that one gets the noticeable result that 
< s z >t >> < s y > T . On the basis of what precedes, Biittiker in refJ 12 ^ intro- 
duced the new Larmor times: r^rp(xi, xf, k) and T^(xi,xf,k), both defined analogously 
to quantities r^f and r^, respectively; as well as the hybrid Larmor times, defined as 
follows: 

(^f R ) 2 = (rilnf + (r 2 L T a , R ) 2 . (18) 

However, the introduction of so many time durations for a single collision (e.g., 
transmission and reflection processes) seems to us physically unjustified. 
Let us notice that, for an opaque rectangular barrier, we obtain 
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(19) 



which results to be different from both the extrapolated Phase Times (12') and the Dwell 
Times (15). 

1-2.6 A complex time approach 

As it is known, a formal generalization of the classical time spent by a particle inside 
the barrier can lead, for E < Vq (when the actual presence of the particle, there, is 
forbidden by classical mechanics), to the introduction of a complex time. More generally, 
an analogous extension of classical time to the quantum domain has been recently proposed 
inl 18 ' (see also, e.g., refs.' 19-21 ', and refs. therein). For one-dimensional motion, following 
refJ 18 !, one of the natural "quantum generalizations" of the classical expression 

r[P(t)] = J** dt dx 5[x - P(t)} (20) 

for the time spent inside a region V — where P(t) is the classical path going from X[(ti) 
to Xf(tf) — is the path-integral average 

r Qu (xi,ti;a; f ,t f ) = < r[P( )] > paths , (21) 

in which P( ) is any arbitrary path between the given end-points. 
For the process relative to Fig.I-1 and to eq.(l), one has:' 5,19 ' 

T Q» = m I ix . T Q" = lh ! dx ^ , (22) 

Jv SQ(x) Jv 5P(x) 

where V is nothing but the interval (xi,X{) [or, in particular, (0, a)]; and 8/5P(x) is a 
functional derivative. In general, quantities are complex; and are connected with 
the Larmor Times by the relations 

R e T T,R = T yT,R i I m T T,R = 1~zT,R ■ (23) 
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Of course, complex time is a useful theoretical tool; even if the ordinary tunnelling- 
times should be real. The physical meaning of the imaginary part is still controversial. t 22 ' 

1-2.7 The Biittiker— Landauer time 

In refsJ 23 J the tunnelling times were studied via a new kind of "gedanken ex- 
perimente", namely by supposing the barrier to possess, besides the ordinary (time- 
independent) part, an additional part oscillating in time: 

V(t) = V + Vi cos cut . (24) 

Since the potential V varies with time, the incident particles — if endowed with electric 
charge (or magnetic moment) — can absorb or emit "modulation quanta" Hcu during the 
tunnelling, which leads to the appearance of sidebands with energies E + nhu; [n = 
±1,±2, ...]. In the first-order approximation in V±, it is enough to consider only the 
neighboring sidebands with energies E ± hu. Biittiker and Landauer did obtain the 
following expressions for the relative sideband intensities 

4"V) = I 1 2 - exp^r-) - if , 

(25) 

where A^ 1 ^ and A$ are the perturbed (sideband) and unperturbed transmission- 
amplitudes, respectively; and similarly for the reflection-amplitudes A^ 1 ^ and A$. In 
eqs.(25) the last equalities (~) hold only for the case of opaque [rectangular] barriers and 
not too high frequencies: i. e., for Tiuu small with respect to both E and Vq — E. Moreover, 
= ma/tiK; r^ L = 2mk/[hK,(K 2 + k 2 )]; and r K = m/tiK 2 . One can see that r^ L is 
identical to the Larmor time r^f as given in eq.(19). 

In our opinion, it is not worthwhile to report about the discussions originated by 
Biittiker-Landauer's approach, since they seem to us as being too technical and insuffi- 
ciently justified; let us only quote, here, the refs.' 24-2 ^ However, two results should be 
mentioned. 

First, Hauge and StovnengW did find a simple connection between, on one side, the 
uj — > limits of A^^^/A^ and A^^Iuj) / A^\ and, on the other side, the complex 
times of eqs.(23): 
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4 ±1} H _ ,Vl r Q,. 4^) _ .^l r Qu m) 

even if the physical meaning of such a connection is not yet very clear. 

Second, it is interesting to recall that Bruinsma and Bak^ 28 ' (see also ref.M) proposed 
the characteristic frequency (t^ 1 ") -1 = TiK/ma to give information about the coupling 
between tunnelling and other accompanying channels, rather than about the intrinsic 
tunnelling times. 



1-3 - ABOUT THE POSSIBILITY OF INTRODUCING CLEAR DEFINI- 
TIONS OF r T AND r R . 

1-3.1 A comment on Hauge and Stovneng's conclusions. 

After having reviewed the main definitions and evaluations of the tunnelling times 
(which we also have presented, and criticized, in Sect. 1-2), the authors of refJ 1 ' concluded 
that no definite, acceptable approach still exists to calculating such tunnelling durations. 
As a necessary but not sufficient condition, to be obeyed by any physically acceptable 
expression of the tunnelling and reflection times tt and Tr, those authors did propose the 
following relation [T = Transmitted, in this case!] 



r Dw = | At |2 Tt + | Ar |2 Tr ? (27) 

which they required to be satisfied by the durations calculated via any method (except 
the dwell one, of course, which ab initio does not separate transmission from reflection 
time). Let us observe that the negative conclusion of ref.W, which is actually the main 
conclusion of that review, is based not only on a criticism of all the previously existing 
approaches (a criticism that we made more complete and even stronger), but also on the 
fact that none of them satisfies condition (27). 

However, relation (27) is unacceptable as a general criterion, since it attributes a 
special role (and meaning) to the Dwell Time r Dw , which on the contrary does not pos- 
sess — in general — the physical meaning of global collision-duration, as we showed in 
Sect.I-2.4. 

In the following Section, we are going to show that it is possible to define (and 
calculate) — in a physically meaningful and self-consistent way — those durations r T 
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and tr. 



1-3.2 A general definition of the collision durations; and Applications to Tun- 
nelling 

A direct, general definition of the collision durations was put forth first by Ohmura' 29 ], 
and then improved — and generalized for finite distances — by us.t 6 ' 8 ] Following refsJ 6 ', the 
transmission and reflection durations < Tt >, < tr > (averaged over the corresponding 
flux densities) can be defined, in the considered case of one-dimensional motion in presence 
of a barrier, as follows [T = Traversal]: 

< r T > = < t(xf) >™ - < t(xi) >f n = 

f^dt JF(x t ,t) JZodt Mx u t) " 

(x t - Xi) < v' 1 > + < Ar T > ; (28) 



J °°dEv\gA T : 


2 i$ h (x u x { , E) 


JS°dEv\gA T \ 


2 



<r R > = <t(x i )>l-<t(x i )> hl = 

J^dttJ^t) jr^dttMxut) 
S^dt Jg(x u t) J^dtMxut) " 

= " — #d^CT s 2,1,1 < " > + < Ate > ' (29) 

which hold when the incoming, reflected and transmitted wave-packets do not interfere: 
i.e., are totally separated in space-time. Quantity g = g(E — E) was defined in eqs.(5), 
(6); while the ordinary Phase Times r^ h , r^ h have been defined in eqs.(10) and (11). 
Moreover, quantities and J™ are defined as in eq.(9), just replacing ip(x, k) of eq.(5) 
by ipn = An e~ tkx and ^ni = ifa = A<? e lkx , respectively. Let us stress that our 
equations (28), (29) do implicitly define also the time delays < At t >, < Atr > 
due to transmission and reflection, respectively; as well as the "average" instants 
< t(xf) >Jj5\ < t(x{) >{[, < t(xi) > in at which the corresponding wave-packets [trans- 
mitted, reflected and initial, respectively] pass through point Xf or X[. 

Notice that for quasi-monochromatic wave packets, i.e., when approximation (8) 
holds, eqs.(28), (29) do directly yield the ordinary Phase-Times r^ h , r^ h , given in 
eqs.(10), (11). 
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However, when x- u X{ are not far from the barrier, then it happens that the in- 
coming, reflected and transmitted wave-packets can interfere. Moreover, the flux density 
J(x, t) does in general change its sign with time; for example, the sign of J(0, t) does 
change from + into — approximately a time % d(arg A-&)/dE after the arrival at x = of 
the initial wave-packet. Therefore, the integrals f^dt t J(x,t) do represent in general 
the algebraic sum of positive and negative quantities, so that the probability densities 

J(x, t) dt 

n.dt j( X ,t) 

are not positive definite and do not possess a direct physical sense. 

Each probability density acquires a physical meaning only during those (partial) 
time-intervals in which the corresponding flux-density J(x, t) does not change its direc- 
tion. As a consequence, the previous integrals are to be split into various integrals, each 
one carried over a partial time-interval such that during it the sign of J(x, t) is only 
positive, or only negative. Afterwards, one will sum over all such contributions. 

In other words, we have to deal only with the positive definite probability densities 

dt J+(x,t) ^ dt J-(x,t) 
f-oo dt J+ (x, t) dt J_ (x, t) 

where J + and J__ represent the positive and negative values of J(x,t), respectively. 

Therefore, we do propose as physically adequate definitions for the average transmis- 
sion time and the average reflection time the following expressions: 



< r T > = < t(x f ) >, - < t(xi) > 



+ 



J^dttJ+fat) JZodtt J + (x u t) 



f^dt J+(x t ,t) 



J^dt J + (x u t) 



(30) 



< r R > = < t(xi) >_- < t(xi) > + = 

IZo dt t J-(Xj,t) _ JZo dtt J+(xj,t) 
J^dt J_(x u t) ~ J^dtJ+ixut) ■ 1 } 

Let us notice that, when Xf > a and X{ — > — oo, equation (30) goes into equation 
(28) since in that case J + (xf,t) = Jr(xf,i) = J(xf,t) and J + (xi,t) = J m {xi,t). Anal- 
ogously, when Xi — > — oo, equation (31) goes into equation (29) since in such a case 
J+(xi,t) = J in (xi,t) and J^(x u t) = J K (x h t). 



14 



What precedes, and in particular eqs.(30), (31), lead us to adopt as suitable, strict 
definitions for the very Tunnelling Time and the Reflection Time at the barrier-front (or 
To-and-Fro Time) the following ones: 



< r tun > = < t(a) >+- < t(0) >+ , 



(32) 



< r to _ fro > = < r R (xi = 0) > = < t(0) > - < t(0) >+ , 



(33) 



where one should recall that the barrier starts at the point x — 0. 

According to us, eqs.(32) and (33) are the correct definitions for the "Tunnelling 
time" < T tU n > and the "Reflection-due-to-the- whole-barrier time" (i.e., the Reflection 
time at the barrier front wall) < r to _f ro >. In conclusion, at variance with the authors 
of review' 1 ^, we think that a positive answer can be given to their question about the 
possibility of a precise, meaningful, univocal definition of the Tunnelling and Reflection 
times; such an answer being provided by our equations (30)-(33). 

Unfortunately, simple analytical expressions for those time-durations in the energy 
representation exist only in particular, limiting cases. In general, even for Gaussian or 
quasi-monochromatic wave packets, calculations can be performed only numerically. Any- 
way, eqs.(30)-(33) can be qualitatively tested in an easy way. 

We are left with the question of the time evolution of wave-packets inside the bar- 
rier: a problem which till now was paid attention to only in ref.^. We shall examine it 
in the coming Section. 

Before going on, let here mention — however — that In the second part of this 
paper we shall show that our definition < tt > for the average transmission time results 
to constitute an improvement with respect to the ordinary dwell-time r Dw formula. 

1-3.3 Time evolution of the tunnelling wave— packets inside the barrier 

In ref.t 7 ' calculations were performed of p(x,t) and J(x,t), at different points x in- 
side the barrier, for a Gaussian wave-packet with an energy spread AE = 0.025 E. 
The results of those calculations are presented in Fig. 1-3, for E = |Vo and na = 5/\/2. 
From it, one can see that the times r p (x), Tj + (x) and tj-(x), taken by the maximum 
of p(x,t), J + (x,t) and \J-(x, t)\, respectively, to penetrate the barrier till the depth 
Ax = x, do not depend linearly on x; and that J(x, t) — inside the barrier — does change 
its sign with time, not very far from the barrier forward wall (0 < x < 0.6 a). 

It is worthwhile to notice that: (i) although the continuity equation dp/dt + 
dJ/dx = goes on holding inside the barrier, nevertheless the equality J = vp 
(which is valid for quasi-monochromatic wave packets outside the barrier) is not valid 
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—not even approximately — inside the barrier; and that: (ii) the effective velocities 
v p = (drp/drr) -1 , vj + = (drj + /dx) _1 and Vj- = (drj_/dx) _1 of the maximum of 
p, J + and | J_|, respectively, not only are non-constant as x varies, but also do not coin- 
cide with each other. 

Passing to the mean velocity of the wave-packet while tunnelling through the whole 
barrier, it can be defined in a natural way as follows: 

^ = (34) 

Let us explicitly notice that, if < r tun > does not increase with a, the "effective" 
speed may become arbitrarily large. This would actually happen when the tunnelling 
time can be expressed by the ordinary Phase Time: < r tun >= r^ h (0, a; E); cf. Sect.I- 
2.3. 

Moreover, we can show that even in general, for rectangular barriers and large val- 
ues of a, quantity < r tun > does not depend practically on a. In fact, it is obviously 

< t(a) >™ = < t(a) > + ; so that the corresponding terms, for Xf = a, become equal in 
eq.(32) and in eqs.(28). Therefore, the same will happen — for quasi- monochromatic wave 
packets — for the corresponding term included in eqs.(12). As a consequence, the term 

< t(a) > + does not depend on a for opaque barriers. As to the second term, < t(0) > + 
, of eq.(32), it differs from < t(0) >i n of eqs.(28) [owing to the effect of interference 
between incoming and reflected waves in the flux J+(x; = 0, t) ] by a quantity that de- 
pends on the reflection time r^ h (0, a; E), on the wave packet time-extension [of the order 
of h/AE], and on the form of the wave-packet. Consequently, also such a "second" term 
does not depend on a for opaque barriers. We can conclude that the Hartmann— 
Fletcher's effect is valid even for our definition (32), so that our approach by 
wave-packets confirms that Q.M. predicts Superluminal tunnelling through 
opaque barriers. 

Let us mention, at this point, that another example of barrier — the inverted os- 
cillator potential — was carefully investigated by Barton, [30] through a slightly different 
formalism. In that paper some interesting new results (which are partially similar to our 
ones) have been put forth in connection with the time evolution of the tunnelling wave- 
packets. In particular, Barton met cases in which, under the barrier, wave-groups with 
lower energy travel faster [cf. also our eq.(28), as well as eq.(12)]. 

It is perhaps worthwhile to add the following observations. The arriving, initial 
wave-packet does interfere with the reflected waves, that start to be generated as soon 
as the packet forward tail reaches the barrier edge; in such a way that (when considering 
the profiles of fluxes J(x, t) before the barrier) the backward tail of J m decreases — for 
destructive interference with Jr — in a larger degree than the forward one. This simu- 
lates an increase of the average speed of the entering-flux profile, J + (x,t). Hence, the 
term < t(x) decreases for negative x ~ 0. In other words, the effective (average) 
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flight-time of the approaching packet from the source to the barrier does decrease. Let 
us now consider what happens inside the barrier, for positive x ~ 0. An analogous in- 
terference effect leads to expect an increase of the (effective) tunnelling time < r tun >; 
which consequently will not coincide with the Phase Time r£ h (0, a; E): not even for 
quasi-monochromatic packets. Finally, it is interesting to note, and easy to recognize, 
that the time-flight decrease before the barrier, and the tunnelling-time increase inside 
the barrier, do exactly cancel each other out, so that the total effect vanishes. In any 
case, let us stress that the fact that the entrance time-instant < t(0) > + is decreased 
by the mentioned "distortion" does not obscure the physical sense of our definition of the 
tunnelling time < r tun >, eq.(32). 

Coming back to the time-evolution of the wave packet inside the barrier, we cannot 
describe it as the quasi-classical motion of a particle. We can visually describe it, on the 
contrary, in terms of the "motions" of the three densities p u (x, t), J+(x,t) and J l }(x, t). 

In particular, we can pictorially say that both profiles [bumps] of the "incoming" 
and the "reflected" flux-densities J+(0,i), J_(0,t) (depicted for x = as a function of 
time: see Fig. 1-3) do repeat themselves — with distortion — for increasing values of the 
penetration depth x. Such a "transmission" takes place rather rapidly (when considering 
the velocities vj + , t>j_ of the bump maxima) in comparison with the initial speed v, 
non- uniformly, and with a gradual decreasing of the bumps [the non-uniformity of such 
"motions" being immediately evident from the non-linear dependence on x of the station- 
ary wave-function phases inside the barrier: cf. eqs.(l)]. In particular, such decreasing is 
so strong for the reflected waves, that they disappear at the end of the barrier. At last, 
we can qualitatively say that the probability density p(x, t) does conserve the form of its 
temporal shape at every position x inside the barrier, but does exponentially decreases as 
the depths increases, and quickly, non-uniformly moves towards the final barrier wall. 



1-4 - CONCLUSIONS OF PART I AND PERSPECTIVES 

In the Fist Part of this work, besides a critical review of the previous definitions of 
the tunnelling times and of their consequences, we proposed for them — on the basis of 
our general formalism' 6,8 ' — new definitions, eqs.(30)-(33), that we regard as physically 
acceptable. We do share with the authors of review' 1 ', however, the opinion that none of 
the previously known definitions were generally acceptable. 

In the First Part above we began analyzing — for the first time — also the time- 
evolution of the tunnelling process inside the barrier. 

Let us recall that our formalism is based on the introduction and application of a 
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[non-selfadjoint, but hermitian] operator for Time as well as on suitable definitions for 
the observables' time-averages [cf. eqs.(28)-(31)]. The physical self-consistency of such a 
formalism may be regarded as supported by results as: (i) the validity of a correspondence 
principle between the time-energy QM commutation relation and the CM Poisson brack- 
ets; (ii) the validity of an Ehrenfest principle for the average time-durations; I 6 ' 8 ' (iii) 
the coincidence of the quasi-classical limit of our own QM definitions for time durations 
(when such a limit exists: i.e., for above-barrier energies) with the analogous, well-known 
expressions of classical mechanics [see in particular the second and third refs.' 6 ', and refs. 
therein]. Moreover, definitions similar to eqs.(28)-(29) have been already applied and 
tested in the analysis of experimental data on nuclear reaction durations, in the range 
1CT 21 — lCT 15 s, obtained by blocking-effect, X-ray spectroscopy and brehmsstrahlung ex- 
periments [see in particular the first two refs.' 6 ', and refs. therein]. Let us stress that 
for completely extracting the time-duration values from experimental data on the above- 
mentioned (non-stationary) processes, which are always non-linearly depending on the 
nuclear reaction durations, it is necessary to have recourse not only to equations of the 
type (28)-(29), but also to correct definitions for the duration variances, and for the dura- 
tion distribution higher-order central moments: ' 6 ' which is provided by our formalism. At 
last, let us mention that such a formalism did provide also useful tools for the resolution 
of some long-standing problems related to the time-energy uncertainty relation. [6, 8] 

As to the future, let us before all recall that analytical expressions in the energy 
representation for < r tun > and < r to _ fro > do not exist, even in the case of the simplest 
barriers. However, numerical calculations can be straightforwardly performed about the 
time evolution of the considered wave-packet inside the barrier. On the basis of such 
one-particle (and one-dimensional) calculations, it will be possible to start developing a 
kinematical theory for the tunnelling of bound or metastable many-particle systems, and 
of unbound aggregates, through various barriers. 

Within the field of nuclear physics, it will be interesting — on the basis of the pe- 
culiarities of the tunnelling-process temporal evolution (Sect. 1-3. 3 above; and ref.' 7 ')— 
to investigate the possibility of observing effects due to the change in form, in volume, 
in orientation and in life-time of many-particle systems (nuclei, fragments,...) during 
their tunnelling. As well as effects due to "collisions", inside the barrier (e.g., the ion- 
ion barrier), between different, successive, penetrating particles: for example, the role of 
those possible collisions in enhancing two-proton sub-barrier transfer in heavy ion inter- 
actions. Another future task will be developing a multi-dimensional description, in terms 
of tunnelling processes, of the sub-barrier fusion of two nuclei: taking account of nuclear 
deformation, formation of noses or a neck, and of dissipation phenomena, f 31 ' 

At last, we should comment on the fact that — when < r tun > does not practically 
depend on the barrier width a, as we have seen to happen for opaque or high barriers 
(Sects. 1-3. 3 and 1-2.3) — then one apparently meets, in connection with eq.(34), speeds 
that (inside the barrier) can assume arbitrarily large values. This does not violate any 
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postulate, as far as we deal with non-relativistic (quantum) physics; and in fact such a 
phenomenon has already been frequently met within "quantum systems". ^ 

It is easy to check that the same would happen, in our case, even when replacing the 
Schroedinger equation by the Klein-Gordon or Dirac equations, at least for barriers that 
do not depend explicitly on time. More interesting is the occurrence of this fact in such 
"quantum field theories" (which deal, however, with semi-classical potentials). 

One could think that those infinitely large speeds might disappear in a self- consistent 
relativistic quantum theoretical treatment. But such an expectation seems to be wrong, 
for reasons that come also from Special Relativity itself. [33] The Hartman (Superluminal 
tunnelling) effect should then be added to the already known results that suggest the 
existence of Superluminal motions: cf., e.g., refs.' 34 '. Results that have been predicted, 
inside relativistic theories like QFT, by many authors, like Sudarshant 35 ! , Van Dam and 
Wigner^ 36 !, Recanti^, Ne'eman' 37 ! and others. 



1-5 - Note added in the 1992 proofs: 

In two recent papers^ 38 ' 39 ! , Biittiker and Landauer published some other critical comments 
about the various approaches examined by us in Sect. 1-2. Let us mention, among them, 
a further criticism of Hauge and Stovneng's relationship [eq.(27), Sect. 1-3.1], appeared in 
ref.' 39 '. Moreover, an interesting, brief review of the first experimental measurements of 
tunnelling times did appear in ref.I 38 '. 



CAPTIONS OF THE FIGURES OF PART I: 

Fig.I-1 - The case of stationary scattering and tunnelling. In this figure we depict a generic 
potential-barrier V(x) (which does not depend explicitly on time), with the three regions 
generated by its presence; and the incoming, reflected and transmitted plane waves. 

Fig. 1-2 - Orientations of spin and magnetic field for the case of the "Larmor clock" (see 
the text). 

Fig. 1-3 - A pictorial view of the time dependence of J(x,t) and p(x,t), for various val- 
ues of the penetration depth x inside a rectangular barrier (for Gaussian wave-packets). 
Actually, as x increases, all the "bumps" suffer an exponential damping; but, for conve- 
nience, we neglected in this figure the exponential factor exp(— kx). 



19 



BIBLIOGRAPHY OF PART I: 



References 

[1] E.H. Hauge and J. A. Stovneng: Rev. Mod. Phys. 61 (1989) 917. 
[2] L. Eisenbud: dissertation (Princeton Univ., 1948). 
[3] D. Bohm: Quantum Theory (Prentice Hall; New York, 1951). 
[4] E.P Wigner: Phys. Rev. 98 (1955) 145. 

[5] F. Capasso, K. Mohammed and A.Y. Cho: IEEE J. Quantum. Electron. 22 (1986) 
1853. 

[6] V.S. Olkhovsky: Sov. J. Part. Nucl. 15 (1984) 130 [English translation]; USSR- 
Doctorate Dissertation (Inst. Nucl. Res., Ukr.Acad.Sc; Kiev, 1986); Nukleonika 35 
(1990) 99-144. An extended summary of the previous works can be found in V.S. 
Olkhovsky: "The study of nuclear reactions by their temporal properties" , in Nuclear 
Reaction Mechanisms, ed by D.Seeliger and H.Kalka (World Scientific; Singapore, 
1991), p.15. 

[7] V.S. Olkhovsky, V.M. Shilov and B.N Zakhariev: Operator Theory: Adv. and Appl. 
46 (1990) 159. [See also A.O. Barut, M. Bozic, Z. Marie and H. Rauch: Z. Phys. 
A238 (1987) 1; A.O. Barut and M. Bozic: Z. Phys. A330 (1988) 319; A335 (1990) 459; 
Physica B151 (1988) 180]. 

[8] See, e.g., V.S. Olkhovsky, E. Recami and A.I. Gerasimchuk: Nuovo Cimento A22 
(1974) 263; E. Recami: "A time operator and the time-energy uncertainty rela- 
tion", in The Uncertainty Principle and Foundations of Quantum Mechanics, ed. by 
W.C.Price and S.S.Chissick (J.Wiley; London, 1977), p. 21; E. Recami: "An opera- 
tor for Time", in Proceedings of the Karpacz Winter School (Recent Developments in 
Relativistic QFT and its Application, vol. 2), ed. by W.Karwowski (Wroclaw Univ. 
Press; Wroclaw, 1976), p. 251, and refs. therein. 

[9] T.E. Hartmann: J. Appl. Phys. 33 (1962) 3427. 

[10] J.R. Fletcher: J. Phys. C18 (1985) L55. 

[11] F.T. Smith: Phys. Rev. 118 (1960) 349. 

[12] M. Biittiker: Phys. Rev. 27 (1983) 6178. 

[13] M.S. Goldberger and K.M. Watson: Collision theory (J.Wiley; New York, 1964). 



20 



[14] J.M. Jauch and J.P. Marchand: Helv. Phys. Acta 40 (1967) 217. 
[15] Ph.A. Martin: Acta Phys. Austriaca Suppl. 23 (1981) 157. 

[16] V.L. Lyuboshitz: Mater. XIX Zymney Shkoly L.I. Ya.F., vol.3 [in Russian] (L.I.Ya.F.; 
Leningrad, 1984), p. 33. 

[17] A.I. Baz': Sov.J.Nucl.Phys. 4 (1967) 182; Sov. J.Nucl.Phys. 5 (1967) 161. 

[18] D. Sokolovski and L.M. Baskin: Phys. Rev. A36 (1987) 4604. 

[19] E. Pollak and W.H. Miller: Phys. Rev. Lett. 53 (1984) 115. 

[20] E. Pollak: J. Chem. Phys. 83 (1985) 1111. 

[21] C.R. Leavens and G.C. Aers: Solid State Commun. 63 (1987) 1101. 

[22] See, e.g., E.Recami, W.A. Rodrigues Jr. and P. Smrz: Hadronic J. 6 (1983) 1773, 
and refs. therein; E.Recami: refs.®; V.S. Olkhovsky and E.Recami: Lett. Nuovo 
Cim. [1st series] 4 (1970) 1165; J. A. Gallardo, A.J. Kalnay and S.H. Risenberg: 
Phys. Rev. 158 (1967) 1484; J.A. Gallardo, B.A. Stec, A.J. Kalnay and B.P Toledo: 
Nuovo Cimento A49 (1967) 393; A48 (1967) 997, 1008. See also H.K. Wimmel: Lett. 
Nuovo Cim. 2 (1971) 363, E674; Nature Phys. Sc. 236 (1972) 79; E. Recami: Riv. 
Nuovo Cim. 9 (1986), issue no. 6, Sect. 13 and refs. therein. 

[23] M. Biittiker and R. Landauer: Phys. Rev. Lett. 49 (1982) 1739; Phys. Scr. 32 (1985) 
429; Ibm J. Res. Dev. 30 (1986) 451. 

[24] S. Collins, D. Lowe and J.R. Barker: J. Phys. C20 (1987) 6213. 

[25] M. Biittiker and R. Landauer: J. Phys. C21 (1988) 6207. 

[26] B. Lowe and S. Collins: J. Phys. C21 (1988) 6210. 

[27] C.R. Leavens: Solid State Commun. 68 (1988) 13. 

[28] R. Bruinsma and P. Bak: Phys. Rev. Lett. 56 (1986) 420. 

[29] T. Ohmura: Prog. Theor. Phys. Suppl. 12 (1964) 108. 

[30] G. Barton: Ann. of Phys. 166 (1986) 322. 

[31] J. Schneider and H. Wolter: private communication (Munich, 1990). 

[32] See, e.g., F. Selleri: Quantum paradoxes, and physical reality (Kluver Ac. Pub.; 
Dordrecht, 1990), and refs. therein; G.C. Hegerfeld: Phys. Rev. Lett. 54 (1985) 
2395; H.P Stapp: Am. J. Phys. 53 (1985) 306. 



21 



[33] See, e.g., E. Recami: Riv. Nuovo Cim. 9 (1986), issue no. 6, pp. 1-178; and refs. 
therein. 

[34] G.C. Hegerfeld: Phys. Rev. Lett. 54 (1985) 2395; E. Recami and W.A. Rodrigues 
Jr.: Prog. Part. Nucl. Phys. 15 (1985) 499, and refs. therein; A.O. Barut, C.K.E. 
Schneider and R. Wilson: J. Math. Phys. 20 (1979) 2244; K. Rafanelli: Phys. Rev. 
D17 (1978) 640; L. Robinett: Phys. Rev. D18 (1978) 3610; A.M. Gleeson, M.G. 
Gundzig, E.C.G. Sudarshan and A. Pagnamenta: Phys. Rev. A6 (1972) 807; ; M. 
Jehle: Phys. Rev. D6 (1972) 441; V.S. Olkhovsky and E.Recami: Nuovo Cimento 
A63 (1969) 814, page 822. 

[35] O.M. Bilaniuk, V.K. Deshpande and E.C.G. Sudarshan: Am. J. Phys. 30 (1962) 618; 
E.C.G. Sudarshan: Ark. Fyz. 39 (1969) 585; J. Dhar and E.C.G. Sudarshan: Phys. 
Rev. 174 (1968) 1808. 

[36] H. Van Dam and E.P Wigner: Phys. Rev. B138 (1965) 1576; 142 (1966) 838. See 
also D.G. Currie, T.F. Jordan and E.C.G. Sudarshan: Rev. Mod. Phys. 35 (1963) 
350. 

[37] Y. Ne'emann: Found. Phys. 16 (1986) 361; in Proc. Symp. on the Foundations of 
Modern Physics, ed. by P. Lahti and P. Mittelstaed (1985), pp. 481; Proc. Nat. 
Acad. Sci. N.Y. 80 (1983) 7051; 77 (1980) 720. 

[38] M. Buttiker: "Traversal, Reflection and Dwell time for quantum tunnelling" , in Elec- 
tronic Properties of Multilayers and Low-dimensional Semiconductors Structures, ed. 
by J. M. Chamberlain et al. (Plenum Press; New York, 1990), pp.297. 

[39] R. Landauer: Ber. Bunsenges. Phys. Chem. 95 (1991) 404. 



22 



SECOND PART 



II- 1. — Introduction 

In the First Part' 11-1 ' of this paper, we have put forth an analysis of the main 
theoretical definitions of the sub-barrier tunnelling and reflection times, and proposed 
new definitions for such durations which seem to be self-consistent within conventional 
quantum mechanics.* 1 In particular, the "prediction" by our theory! 11-1 ^ of the reality 
of the Hariman effect^ n ~ 2 ^ in tunnelling processes has recently received — due to the 
analogy' 11-3 ] between tunnelling electrons and evanescent waves — quite interesting, even 
if indirect, experimental verifications at Cologne, t 11- ^ Berkeley, P 1-5 ' Florence ' II_6 ' and 
Vienna. [II " 6] 

Main aims of this Second Part [that appeared in print in J. de Physique-I 5 (1995) 
1351] are: presenting and analysing the results of several numerical calculations of the 
penetration and return times inside a rectangular potential barrier during tunnelling 
(Sect.II-3); and proposing new suitable formulae for the distribution variances of the 
transmission and reflection times (Sect.II-2). The results of our numerical evalua- 
tions seem to confirm that our approach is physically acceptable, and that it moreover 
implied, and implies, the existence of the so-called "Hariman effect" (i.e., of tunnelling 
with Superluminal group-velocities) even for non-quasi-monochromatic packets. 

Before all, let us add here — however — some brief comments about a few furthet 
papers, appeared recently: 

(i) First, let us mention that in Part I above we have overlooked a new expression 
for the dwell-time r Dw derived by Jaworsky and WardlawP 1-7 ' 8 ' 

r^(x u x ( ;k) = (f^dttJixut) - dtt J(x u t)) dt J ta (x u t))' 1 , (II-l) 

which is indeed equivalent t 11-7 ' to our eq.(16) of ref.[II-l] (all notations being defined 
therein): 

T^(x u xf,k) = (jZodtJ^dxpix^^dtJUxut))' 1 . (II-2) 

This equivalence reduces the difference, between our definition < r T > of the av- 
erage transmission time — under our assumptions — and quantity r Dw , to the dif- 

0#1 Let us take advantage of the present opportunity for pointing out that a misprint entered our 
eq.(10) in ref.[II-l], whose last term ka ought to be eliminated. Moreover, due to an editorial error, in 
the footnote at page 32 of our rcf. [11-16] the dependence of G on Afc disappeared, whilst in that paper 
we had assumed G(k - k) = C exp[-(fc - fc) 2 /(Afc) 2 ]. 
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ference between the average made by using the positive-definite probability density 
dt J+(x, t) I /f^ dt J+(x, t) and the average made by using the ordinary "probability den- 
sity" dtJ(x,t)/ Jf^ dt J(x, t). Generally speaking, the last expression is not always 
positive definite, as it was explained at page 350 of ref.[II-l], and hence does not possess 
any direct physical meaning. 

(ii) In ref. [II- 9] an attempt was made to analyze the evolution of the wave packet 
mean position < x(t) > ("center of gravity"), averaged over pdx, during its tunnelling 
through a potential barrier. Let us here observe that the conclusion to be found therein, 
about the absence of a causal relation between the incident space centroid and its trans- 
mitted equivalent, holds only when the contribution coming from the barrier region to the 
space integral is negligible. 

(iii) Let us also add that in ref. [11-10] it was analyzed the distribution of the trans- 
mission time tt in a rather sophisticated way, which is very similar to the dwell-time 
approach, however with an artificial, abrupt switching on of the initial wave packet. We 
are going to propose, on the contrary, and in analogy with our eqs.(30)-(31) in ref. [II- 1], 
the following expressions, as physically adequate definitions for the variances (or disper- 
sions) D t t and Dr R of the transmission and reflection time [see Sect.II-2], respectively: 

Dt t = Dt+(x t ) +Dt + (xi) (11-3) 

and 

Dt r e D*_(xi)+Dt_(a;i) , (II-4) 

where 

JZodt 1?J ± (x,t) ( JZodttJ ± (x,t) 
±{X> ~ SZodtJ ± (x,t) [ JZodtJ^t) 

Equations (1-3)— (1-5) are based on the formalism expounded in refs. [11-11], as well as on 
our definitions for J±(x,t) in ref.[II-l]. Of course, we are supposing that the integrations 
over J + (xf)dt, J + (x[) dt and J^(xi)dt are independent of one another. We shall 
devote Sect.II-2, below, to these problems, i.e., to the problem of suitably defining mean 
values and variances of durations, for various transmission and reflection processes during 
tunnelling. 

(iv) Below, in Sect.II-4, we shall briefly re-analyse some other definitions of tun- 
nelling durations. 

Before going on, let us recall that several reasons "justify" the existence of different 



(II-5) 
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approaches to the definition of tunnelling times: (a) the problem of defining tunnelling 
durations is closely connected with that of defining a time operator, i.e., of introducing 
time as a (non-selfadjoint) quantum mechanical observable, and subsequently of adopt- 
ing a general definition for collision durations in quantum mechanics. Such preliminary 
problems did receive some clarification in recent times (see, for example, ref.[II-l] and 
citations [8] and [22] therein); (b) the motion of a particle tunnelling inside a potential 
barrier is a purely quantum phenomenon, devoid of any classical, intuitive limit; (c) the 
various theoretical approaches may differ in the choice of the boundary conditions or in 
the modelling of the experimental situations. 

II-2. — Mean values and Variances for various Penetration and Return Times 
during tunnelling 

In our previous papers, and in Part I above, we proposed for the transmission and 
reflection times some formulae which imply — as functions of the penetration depth- 
integrations over time of J+(x, t) and J_(x,t), respectively. Let us recall that the total 
flux J(x, t) inside a barrier consists of two components, J + and J_, associated with motion 
along the positive and the negative x-direction, respectively. Work in similar directions 
did recently appear in ref. [11-12] . 

Let us refer ourselves — here — to tunnelling and reflection processes of a particle 
by a potential barrier, confining ourselves to one space dimension. Namely, let us study 
the evolution of a wave packet ^(x, t), starting from the initial state ty m (x,t); and follow 
the notation introduced in ref. [II- 1]. In the case of uni-directional motions it is already 
knownt 11 " 13 ! that the flux density J(x,t) = Re[(ih/m)^(x,t) d^*(x,t)/dx] can be 
actually interpreted as the probability that the particle (wave packet) passes through 
position x during a unitary time-interval centered at t, as it easily follows from the 
continuity equation and from the fact that quantity p(x,t) = \^(x, t)\ 2 is the 
probability density for our "particle" to be located, at time t, inside a unitary space- 
interval centered at x. Thus, in order to determine the mean instant at which a moving 
wave packet ty(x,t) passes through position x, we have to take the average of the time 
variable t with respect to the weight w(x, t) = J(x, t)/ J(x, t) dt. 

However, if the motion direction can vary, then quantity w(x, t) is no longer positive 
definite, and moreover is not bounded, because of the variability of the J(x, t) sign. In 
such a case, one can introduce the two weights: 




(II - 6) 




(II - 7) 



25 



where J + (x,t) and J-(x,t) represent the positive and negative parts of J(x,t), respec- 
tively, which are bounded, positive-definite functions, normalized to 1. Let us show 
that, from the ordinary probabilistic interpretation of p(x,t) and from the well-known 
continuity equation 

dp(x,t) dJ(x,t) 

~^ + ^x—-^ (II " 8) 
it follows also in this (more general) case that quantities w + and W-, represented by 
eqs.(II-6), (H-7), can be regarded as the probabilities that our "particle" passes through 
position x during a unit time-interval centered at t (in the case of forward and backward 
motion, respectively). 

Actually, for those time intervals for which J = J + or J = J_, one can rewrite 
eq.(II-8) as follows: 

dp>(x,t) dJ+(x,t) 

~^r~ ~ dx~ (n " 9M) 

dp^Xit) dJ-(x,t) 



dt dx 



(II -9.6) 



respectively. Relations (9. a) and (9.b) can be considered as formal definitions of dp>/dt 
and dp</dt. Let us now integrate eqs.(II-9.a), (II-9.b) over time from — oo to t; we 
obtain: 



,t) = -f t d J±^ll d t> (II-10.a) 

J-oo OX 



p>(x 



P< (x,t) = -f t °±=pl±dl> (II -10.6) 

J-oo OX 

with the initial conditions p>(x, — oo) = p<(x, — oo) = 0. Then, let us introduce the 
quantities 

/•OO rt 

A^>(a;,oo;t) = / p>(x' ,t) dx' = / J+(x,t')dt' >0 (II - ll.a) 



(-oo,x;t)= f P< {x',t)dx' = - f J_(x,t')dt' >0, (II -11.6) 

J —oo J —oo 
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which have the meaning of probabilities for our "particle" to be located at time t 
on the semi-axis (x, oo) or (— oo,x) respectively, as functions of the flux densities 
J + (x,t) or J_(x,t), provided that the normalization condition /f^ p(x, t)dx = 1 is 
fulfilled. The r.h.s.'s of eqs.(II-ll.a) and (11-11. b) have been obtained by integrat- 
ing the r.h.s.'s of eqs.(II-lO.a) and (11-10. b) and adopting the boundary conditions 
J+(— oo,t) = J-(— oo,t) = 0. Now, by differentiating eqs.(II-ll.a) and (11-11. b) with 
respect to t, one obtains: 

d^ix.oo.t) , , 

— = J+(x,t) > (II-12.a) 



dt 

dN<(x, — oo, t) 



-J-(x,t)>0. (II -12.6) 



dt 

Finally, from eqs.(II-ll.a), (11-11. b), (11-12. a) and (11-12. b), one can infer that: 

dN>(x, oo; t)/dt 

w + (x, t) = — —. r- (II -13. a) 

v ' N>(x, -oo; oo) v ' 

, . dN < (x,—oo;t)/dt . ,„,. 

w-(x, t) = J) , II - 13.6) 

v ' iV<(-oo,:r; oo) v ' 

which justify the abovementioned probabilistic interpretation of w+(x,t) and w-(x,t). 
Let us notice, incidentally, that our approach does not assume any ad hoc postulate, con- 
trarily to what believed by the author of ref. [11-14]. 

At this point, we can eventually define the mean value of the time at which our 
"particle" passes through position x, travelling in the positive or negative direction of the 
x axis, respectively, as: 

<i+W>= WMS (II " 14a) 

and, moreover, the variances of the distributions of these times as: 

Dt+W " -■<*»<*»'' < n - 15 - fl > 

Di - W " ¥Wm - [<t - (x)>f > (II - 15 - 6) 
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in accordance with the proposal presented in refs. [11-1,15]. 

Thus, we have a formalism for defining mean values, variances (and other central 
moments) related to the duration distributions of all possible processes for a particle, 
tunnelling through a potential barrier located in the interval (0, a) along the x axis; and 
not only for tunnelling, but also for all possible kinds of collisions, with arbitrary energies 
and potentials. For instance, we have that 

< T T (x i7 Xf) > = < t+(Xf) > - < t+(Xi) > (II - 16) 

with — oo < Xi < and a < Xf < oo; and therefore (as anticipated in eq.(II-3)) that 

DT T (x u x f ) = Dt + (xf) +Dt + (xi) , 

for transmissions from region (— oo, 0) to region (a, oo) which we called t 11 " 1 ' regions I and 
III, respectively. Analogously, for the pure (complete) tunnelling process one has: 

< r Tun (0, a) > = < t + (a) > - < t+(0) > (II - 17) 

and 

Dr Tun (0,a) ee Dt + (a) + Dt + (0) ; (11-18) 

while one has 

< r Pen (0,x f ) > ee < t+(x f ) > - < i+(0) > 

and 

Dr Pcn (0,x f ) ee Dt + (x f ) + Dt + (0) (11-20) 

[with < Xf < a] for penetration inside the barrier region (which we called region II). 
Moreover: 

< T Ret (x,x) > ee < t-(x) > - < t + (x) > (II -21) 

Dr Rct (x,x) ee Dt4x) + Dt + (x) (11-22) 

[with < x < a] for ''return processes' 11 inside the barrier. At last, for reflections in 
region I, we have that: 

< T R (X U Xi) > EE < t-(Xi) > - < t+(Xi) > (II - 23) 



(11-19) 
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[with — oo < X{ < a], and (as anticipated in eq.(II-4)) that DT^(xi,Xi) = Dt_(x;) + 
Dt+(xi). 

Let us stress that our definitions hold within the framework of conventional quantum 
mechanics, without the introduction of any new postulates, and with the single measure 
expressed by weights (13. a), (13.6) for all time averages (both in the initial and in the 
final conditions). 

According to our definition, the tunnelling phase time (or, rather, the transmission 
duration), defined by the stationary phase approximation for quasi- monochromatic parti- 
cles, is meaningful only in the limit x\ — > oo when J + (x, t) is the flux density of the initial 
packet Ji n of incoming waves (in absence of reflected waves, therefore). 

Analogously, the dwell time, which can be represented (cf. eqs.(II-l),(II-2)) by the 



expression 



[H-7,8,16] 



T Dw (Xi,X f ) 



* J(x f ,t) dt- t J(x u t) dt 



/oo 
-oo 



dt 



with — oo < X[ < 0, and Xf > a, is not acceptable, generally speaking. In fact, the weight 
in the time averages is meaningful, positive definite and normalized to 1 only in the rare 
cases when x\ — > — oo and J m = Jm (i.e., when the barrier is transparent). 



II-3. — Penetration and Return process durations, inside a rectangular barrier, 
for tunnelling gaussian wave packets: Numerical results 

We put forth here the results of our calculations of mean durations for various pen- 
etration (and return) processes, inside a rectangular barrier, for tunnelling gaussian wave 
packets; one of our aims being to investigate the tunnelling speeds. In our calculations, 
the initial wave packet is 

roo 

*in(z, t) = G(k- k) exp[iA;x - iEt/h] dk (II - 24) 

Jo 

with 

G(k-k) =Cexp[-(k-k) 2 /(2Ak) 2 ] , (11-25) 

exactly as in ref. [II-8]; and with E = h 2 k 2 /2m; quantity C being the normalization 
constant, and m the electron mass. Our procedure of integration was described in ref. [II- 
16]. 

Let us express the penetration depth in angstroms, and the penetration time in 
seconds. In Fig.II-1 we show the plots corresponding to a = 5 A, for AA; = 0.02 and 
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0.01 A , respectively. The penetration time < r Pcn > always tends to a saturation 
value. 

In Fig.II-2 we show, for the case Ak = 0.01 A , the plot corresponding to a = 
10 A. It is interesting that < rp cn > is practically the same (for the same Ak) for a = 5 
and a = 10 A, a result that confirm, let us repeat, the existence^ 11 " 1 ] of the so-called 
Hartman effect} 11 ' 2 ^ Let us add that, when varying the parameter Ak between 0.005 
and 0.15 A and letting a to assume values even larger than 10 A, analogous results 
have been always gotten. Similar calculations have been performed (with quite reasonable 
results) also for various energies E in the range 1 to 10 eV. #2 

In Figs.II-3, 4 and 5 we show the behaviour of the mean penetration and return 
durations as function of the penetration depth (with x- x = and < X{ = x < a), for 
barriers with height Vq = 10 eV and width a = 5 A or 10 A. In Fig.II-3 we present the 
plots of < rp en (0,x) > corresponding to different values of the mean kinetic energy: E 
= 2.5 eV, 5 eV and 7.5 eV (plots 1, 2 and 3, respectively) with Ak = 0.02A; and E = 5 
eV with Ak = 0.04A (plot 4), always with a = 5A. In Fig.II-4 we show the plots 
of < r Pcn (0,:r) >, corresponding to a = 5A, with Ak = 0.02A and 0.04 A (plots 1 
and 2, respectively); and to a = 10A, with Ak = 0.02A" 1 and 0.04 A" 1 (plots 3 and 4, 
respectively), the mean kinetic energy E being 5 eV, i.e., one half of Vq. In Fig.II-5 the 
plots are shown of < T-& e t(x, x) >. The curves 1, 2 and 3 correspond to E — 2.5 eV, 5 eV 
and 7.5 eV, respectively, for Ak = 0.02A and a = 5A; the curves 4, 5 and 6 correspond 
to E = 2.5 eV, 5 eV and 7.5 eV, respectively, for Ak = 0.04A and a = 5A; while the 
curves 7, 8 and 9 correspond to Ak = 0.02A and 0.04 A , respectively, for E = 5 eV 
and a = lOA. 

Also from the new Figs.II-3-5 one can see that: 1) at variance with ref.[II-8], no 
plot considered by us for the mean penetration duration < r Pcn (0,x) > of our wave 
packets presents any interval with negative values, nor with a decreasing < rp en (0,x) > 
for increasing x; and, moreover, that 2) the mean tunnelling duration < TT un (0,a) > 
does not depend on the barrier width a ("Hartman effect"); and finally that 3) quantity 
< TTun(0,a) > decreases when the energy increases. Furthermore, it is noticeable 
that also from Figs.II-3-5 we observe: 4) a rapid increase for the value of the electron 
penetration time in the initial part of the barrier region (near x — 0); and 5) a tendency 
of < Tp en (0,x) > to a saturation value in the final part of the barrier, near x — a. 

o#2 p or interested reader, let us recall that, when integrating over dt, we used the interval 
— 10~ 13 s to +10~ 13 s (symmetrical with respect to t = 0), very much larger than the tempo- 
ral wave packet extension. [Recall that the extension in time of a wave packets is of the order of 

l/(vAk) = (Ak ^JiE/m)- 1 ~ 1CT 16 s]. Our "centroid" has been always t = 0; x a = 0. For 
clarity's sake, let us underline again that in our approach the initial wave packet ^m^, t) is not regarded 
as prepared at a certain instant of time, but it is expected to flow through any (initial) point Xi during 
the infinite time interval (— oo, +oo), even if with a finite time-centroid to. The value of such centroid 
to is essentially defined by the phase of the weight amplitude G(k — k), and in our case is equal to when 
G(k — k) is real. 



30 



Feature 2), firstly observed for quasi-monochromatic particles, t 11-2 ' does evidently 
agree with the predictions made in ref.[II-l] for arbitrary wave packets. Feature 3) is also 
in agreement with previous evaluations performed for quasi-monochromatic particles and 
presented, for instance, in refs. [11-1,2,15]. Features 4) and 5) can be apparently explained 
by interference between those initial penetrating and returning waves inside the barrier, 
whose superposition yields the resulting fluxes J+ and J_. In particular, if in the initial 
part of the barrier the returning- wave packet is comparatively large, it does essentially 
extinguish the leading edge of the incoming- wave packet. By contrast, if for growing x 
the returning-wave packet quickly vanishes, then the contribution of the leading edge of 
the incoming-wave packet to the mean penetration duration < rp cn (0,x) > does initially 
(quickly) grow, while in the final barrier region its increase does rapidly slow down. 

Furthermore, the larger is the barrier width a, the larger is the part of the back edge 
of the incoming-wave packet which is extinguished by interference with the returning- 
wave packet. Quantitatively, these phenomena will be studied elsewhere. Finally, in 
connection with the plots of < r Rct (:r, x) > as a function of x, presented in Fig.II-5, let 
us observe that: (i) the mean reflection duration < r R (0,0) > = < r Ret (0, 0) > does 
not depend on the barrier width a; (ii) in correspondence with the barrier region betwen 
and approximately 0.6 a, the value of < r Rct (0,a;) > is almost constant; while (iii) 
its value increases with x only in the barrier region near x = a (even if it should be 
pointed out that our calculations near x = a are not so good, due to the very small 
values assumed by ff^ J-{x, t)dt therein). Let us notice that point (i), also observed 
firstly for quasi-monochromatic particles, P 1 " 2 ' is as well in accordance with the results 
obtained in ref.[II-l] for arbitrary wave packets. Moreover, also points (ii) and (iii) can 
be explained by interference phenomena inside the barrier: if, near x — a, the initial 
returning-wave packet is almost totally quenched by the initial incoming-wave packet, 
then only a negligibly small piece of its back edge (consisting of the components with the 
smallest velocities) does remain. With decreasing x (x — > 0), the unquenched part of 
the returning-wave packet seems to become more and more large (containing more and 
more rapid components), thus making the difference < T^ et (0,x) > — < Tp en (0,x) > 
almost constant. And the interference between incoming and reflected waves at points 
x < does effectively constitute a retarding phenomenon [so that < t_(x = 0) > is 
larger than < tr(x = 0)], which can explain the larger values of < tr(x = 0, x — 0) > 
in comparison with < T-r un (x — 0, x — a) >. 

Therefore our evaluations, in all the cases considered above, appear to confirm our 
previous analysis at page 352 of ref. [II-l], and our conclusions therein concerning in par- 
ticular the validity of the Hartman effect also for non-quasi-monochromatic wave packets. 
Even more, since the interference between incoming and reflected waves before the barrier 
(or between penetrating and returning waves, inside the barrier, near the entrance wall) 
does just increase the tunnelling time as well as the transmission times, we can expect 
that our non-relativistic formulae for < r Tun (0,a) > and < Tt(x{ < 0,X{ > a) > will 
always forward positive values.* 3 
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At this point, it is necessary — however — to observe the following. Even if our 
non-relativistic equations are not expected (as we have just seen) to yield negative times, 
nevertheless one ought to bear in mind that (whenever it is met an object, O, 
travelling at Superluminal speed) negative contributions should be expected 
to the tunnelling times: and this ought not to be regarded as unphysical. In 
fact, whenever an "object" O overcomes the infinite speech 11-18 ] with respect to a certain 
observer, it will afterwards appear to the same observer as its "anti-object" O travelling 
in the opposite space direction' 11 " 18 ]. For instance, when passing from the lab to a frame 
T moving in the same direction as the particles or waves entering the barrier region, the 
objects O penetrating through the final part of the barrier (with almost infinite speeds, 
like in Figs.II-1-5) will appear in the frame T as anti-objects O crossing that portion 
of the barrier in the opposite space-direction^ 1 ^. In the new frame J 7 , therefore, such 
anti-objects O would yield a negative contribution to the tunnelling time: which could 
even result, in total, to be negative. For any clarifications, see refs. [11-18]. So, we have 
no objections a priori against the fact that Leavens can find, in certain cases, negative 
values' 11 " 8 ' 1 ^: e.g, when applying our formulae to wave packets with suitable initial con- 
ditions. What we want to stress here is that the appearance of negative' 11-2 ' times (it 
being predicted by Relativity itself, I 11-18 ! when in presence of anything travelling faster 
than c) is not a valid reason to rule out a theoretical approach. 

At last, let us — incidentally — recall and mention the following fact. Some pre- 
liminary calculations of penetration times (inside a rectangular barrier) for tunnelling 
gaussian wave packets had been presented by us in 1994 in ref. [11-16]. Later on — looking 
for any possible explanations for the disagreement between the results in ref. [II-8] and in 
our ref. [11-16] — we discovered, however, that an exponential factor was missing in a term 
of one of the fundamental formulae on which the numerical computations (performed by 
our group in Kiev) were based: a mistake that could not be detected, of course, by our 
careful checks about the computing process. Therefore, the new results of ours appearing 
here in Figs.II-1-2 [and appeared in J. de Physique-I 5 (1995) 1351] should replace Figs.l- 
3 of ref. [11-16]. One may observe that, by using the same parameters as (or parameters 
very near to) the ones adopted by for the Figs. 3 and 4 of ref. [II-8], our new, corrected 
figures II-l and 2 result to be more similar to Leavens' than the uncorrected ones (and 

°# 3 A different claim by Delgado, Brouard and Muga^ 11-17 ! does not seem to be relevant to our 
calculations, since it is based once more, like ref. [II-8], not on our but on different wave packets (and 
over-barier components are also retained in ref. [11-17], at variance with us). Moreover, in their classical 
example, they overlook the fact that the mean entrance time < t+(0) > gets contribution mainly by 
the rapid components of the wave packet; they forget, in fact, that the slow components are (almost) 
totally reflected by the initial wall, causing a quantum-mechanical reshaping that contributes to the 
initial "time decrease" discussed by us already in the last few paragraphs of page 352 in ref. [II- 1]. All 
such phenomena reduce the value of < t + (0) >, and we expect it to be (in our non-relativistic treatment) 
less than < t+(a) >. 
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this is of course a welcome step towards the solution of the problem). One can verify- 
once more, however, that our theory appears to yield for those parameters non-negative 
results for < rp en (xf) >, contrarily to a claim in ref.[II-8]. Actually, our previous general 
conclusions have not been apparently affected by the mentioned mistake. In particular, 
the value of < r Pcn (:rf) > increases with increasing Xf, and tends to saturation for Xf — > a. 
We acknowledge, however, that the difference in the adopted integration ranges (—00 
to +00 for us, and to +00 for ref.[II-8]) does not play an important role, contrarily 
to our previous belief, I 11-16 ! in explaining the remaining discrepancy between our results 
and those in ref.[II-8]. Such a discrepancy might perhaps depend on the fact that the 
functions to be integrated do fluctuate heavily* 4 (anyway, we did carefully check that 
our own elementary integration step in the integration over dk was small enough in order 
to guarantee the stability of the numerical result, and, in particular, of their sign, for 
strongly oscillating functions in the integrand). More probably, the persisting disagree- 
ment can be merely due to the fact — as recently claimed also by Delgado et al.' 11-17 '— 
that different initial conditions for the wave packets were actually chosen in ref. [II-8] and 
in ref. [II- 1]. Anyway, our approach seems to get support, at least in some particular 
cases, also by a recent article by Brouard et al., which "generalizes" — even if starting 
from a totally different point of view — some of our results J 11-12,15 '. 

Let us take advantage of the present opportunity for answering other criticisms ap- 
peared in ref. [H-8] , where it has been furthermore commented about our way of performing 
actual averages over the physical time. We cannot agree with those comments: let us re- 
emphasize in fact that, within conventional quantum mechanics, the time t(x) at which 
our particle (wave packet) passes through the position x is "statistically distributed" 
with the probability densities dtJ±(x,t)/ J^° 00 dtJ±(x,t), as we explained at page 350 
of ref. [II- 1]. This distribution meets the requirements of the time-energy uncertainty 
relation. 

We also answered in Sect. II- 1 the comments in ref.[II-8] about our analysis^ 11-1 ' of 
the dwell-time approaches. I 11-19 ! 

The last object of the criticism in ref.[II-8] refers to the impossibility, in our approach, 
of distinguishing between "to be transmitted" and "to be reflected" wave packets at 
the leading edge of the barrier. Actually, we do distinguish between them; only, we 
cannot — of course — separate them, due to the obvious superposition (and interference) 
of both wave functions in p(x,t), in J(x,t) and even in J±(x,t). This is known to be 
an unavoidable consequence of the superposition principle, valid for wave functions in 
conventional quantum mechanics. That last objection, therefore, should be addressed 
to quantum mechanics, rather then to us. Nevertheless, Leavens' aim of comparing 
the definitions proposed by us for the tunnelling times not only with conventional, but 
also with non-standard quantum mechanics might be regarded a priori as stimulating and 

°# 4 We can only say that we succeeded in reproducing results of the type put forth in ref.[II-8] by 
using larger steps; whilst the "non-causal" results disappeared — in the considered cases — when adopting 
small enough integration steps. 
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possibly worth of further investigation. 

II-4. — Further remarks 

In connection with the question of "causality" for relativistic tunnelling particles, 
let us stress that the Hartman-Fletcher phenomenon (very small tunnelling durations), 
with the consequence of Superluminal velocities for sufficiently wide barriers, was found 
theoretically also in QFT for Klein-Gordon and Dirac equations, t 11 " 1 ' and experimentally 
for electromagnetic evanescent-mode wave packets' 11 " 4-6 ! (tunnelling photons). It should 
be recalled that the problem of Superluminal velocities for electromagnetic wave packets in 
media with anomalous dispersion, with absorption, or behaving as a barrier for photons 
(such as regions with frustrated internal reflection) has been present in the scientific 
literature since long (see, for instance, quotations [2,1,18], and refs. therein); even if a 
complete settlement of the causal problem — already available for point particles' 11 " 18 '— 
does not seem to be yet available for the case of relativistic waves. Apparently, it is 
not sufficient to pay attention only to group velocity and mean duration for a "particle" 
passing through a medium; on the contrary, it is important taking into account and 
studying ab initio the variances (and the higher order central moments) of the duration 
distributions, as well as the wave packet reshaping in presence of a barrier, or inside 
anomalous media (even if reshaping does not play always an essential role). 

Passing to the approaches alternative to the direct description of tunnelling processes 
in terms of wave packets, let us here recall those ones which are based on averaging over the 
set of all dynamical paths (through the Feynman path integral formulation, the Wigner 
distribution method, and the non-conventional Bohm approach), and others that use 
additional degrees of freedom which can be used as "clocks" . General analyses of all such 
alternative approaches can be found in refs. [11-1,20-24] from different points of view. 

If one confines himself within the framework of conventional quantum mechanics, 
then the Feynman path integral formulation seems to be adequate J 11 " 24 ' But it is not 
clear what procedure is needed to calculate physical quantities within the Feynman-type 
approach I 11-23 !, and usually such calculations result in complex tunnelling durations. The 
Feynman approach seems to need further modifications if one wants to apply it to the time 
analysis of tunnelling processes, and its results obtaind up to now cannot be considered 
as final. 

As to the approaches based on introducing additional degrees of freedom as "clocks" , 
one can often realize that the tunnelling time happens to be noticeably distorted by the 
presence of such degrees of freedom. For example, the Biittiker-Landauer time is con- 
nected with absorption or emission of modulation quanta (caused by the time-dependent 
oscillating part of the barrier potential) during tunnelling, rather than with the tun- 
nelling process itself. I 11-1 ' 15 ! And, with reference to the Larmor precession time, it has 
been shown' 11 " 11 ' 20 ' that this time definition is connected not only with the intrinsic tun- 
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nelling process, but also with the geometric boundaries of the magnetic field introduced 
as a part of the clock: for instance, if the magnetic field region is infinite, one ends up 
with the phase tunnelling time, after an average over the (small) energy spread of the 
wave packet. Actually, those "clock" approaches, when applied to tunnelling wave pack- 
ets, seem to lead — after eliminating the distortion caused by the additional degrees of 
freedom — to the same results as the direct wave packet approach, whatever be the weight 
function adopted in the time integration. 

A more pedagogical, and detailed, review on the same subject (in Italian, however) 
has been "published" electronically, as LANL Archive ffi |cond-mat/9802126 . 

Before closing this paper, in which we met Superluminal motions, we would like to 
put forth and comment the following information. Since the pioneering work by Bateman, 
it is known that the relativistic wave equations — scalar, electromagnetic and spinor — 
admit solutions with subluminal (v < c) group velocities [11-25]. More recently, also 
Superluminal (v > c) solutions have been constructed for those homogeneous wave equa- 
tions, in refs. [11-26] and quite independently in refs. [11-27]: in some cases just by applying 
a Superluminal Lorentz transformation[II-18,28]. Exactly the same happens in the case 
of acoustic waves, with the presence of "sub-sonic" and "Super-sonic" solutions [11-29]; 
so that one can expect they to exist, e.g., also for seismic wave equations. Or, rather, 
we might expect the same to be true even in the case of gravitational waves. At last, 
it is interesting to remark that some, at least, of the Super-sonic (and Super-luminal) 
solutions, when experimentally realized [11-30], appear to be X-shaped, so as predicted in 
1982 in ref. [11-31]. 

A brief review on the experimental data, that — in four different sectors of physics — 
seem to indicate the existence of Superluminal motions, appeared as an Appendix to the 
paper in LANL Archives # |physics/9712051| (to be published in Physica A, 1998). 
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Captions of the Figures of Part II 

Fig.II-1 - Behaviour of the average penetration time < rp en (0,x) > (expressed in sec- 
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onds) as a function of the penetration depth Xf = x (expressed in angstroms) through 
a rectangular barrier with width a = 5 A, for AA; = 0.02 A (dashed line) and 
AA; = 0.01 A (continuous line), respectively. The other parameters are listed in foot- 
note #1. It is worthwhile to notice that < rp en > rapidly increases for the first, few initial 
angstroms (~ 2.5 A), tending afterwards to a saturation value. This seems to confirm the 
existence of the so-called "Hartman effect" J 11-2 ' 1,15 ] 

Fig.II-2 - The same plot as in Fig.II-1, for AA; = 0.01 A , except that now the barrier 
width is a = 10 A. Let us observe that the numerical values of the (total) tunnelling time 
< tt > practically does not change when passing from a = 5 A to a = 10 A, again in 
agreement with the characteristic features 1 of the Hartman effect. Figures II-l and 2 do 
improve (and correct) the corresponding ones, preliminarily presented by us in ref. [11-16]. 

Fig.II-3 - Behaviour of < Tp en (0, x) > (expressed in seconds) as a function of x (expressed 
in angstroms), for tunnelling through a rectangular barrier with width a = 5 A and for 
different values of E and of AA;: 

curve 1: AA; = 0.02A -1 and E = 2.5 eV; curve 2: AA; = 0.02A" 1 and E = 5.0 eV; 
curve 3: AA; = 0.02A" 1 and E = 7.5 eV; curve 4: AA; = 0.04A" 1 and E = 5.0 eV. 

Fig.II-4 - Behaviour of < r Pcn (0,a;) > (in seconds) as a function of x (in angstroms) for 
E = 5 eV and different values of a and A A;: 

curve 1: a = 5 A and AA; = 0.02A ; curve 2: a = 5 A and AA; = 0.04A ; curve 3: 
a = 10 A and AA; = 0.02A" 1 ; curve 4: a = 10 A and AA; = 0.04A" 1 . 

Fig.II-5 -Behaviour of < r Rct (x, x) > (in seconds) as a function of x (in angstroms) for 
different values of a, E and AA;: 

curve 1: a = 5 A, E = 2.5 eV and AA; = 0.02A" 1 ; curve 2: a = 5 A, E = 5.0 eV and 
AA; = 0.02A" 1 ; curve 3: a = 5 A, E = 7.5 eV and AA; = 0.02A" 1 ; curve 4: a = 5 A, 
E = 2.5 eV and AA; = 0.02A" 1 ; curve 5: a = 5 A, E = 5.0 eV and AA; = 0.02A" 1 ; 
curve 6: a = 5 A, E = 7.5 eV and AA; = 0.02A ; curve 7: a = 5 A, E — 5.0 eV and 
AA; = 0.02A" 1 ; curve 8: a = 5 A, E = 5.0 eV and AA; = 0.02A~\ 
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